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We study the carrier capture and relaxation due to Coulomb scattering in group-Ill nitride quan- 
tum dots on the basis of population kinetics. For the states involved in the scattering processes the 
combined influence of the quantum-confined Stark effect and many-body renormalizations is taken 
into account. The charge separation induced by the built-in field has important consequences on 
the capture and relaxation rates. It is shown that its main effect comes through the renormalization 
of the energies of the states involved in the collisions, and leads to an increase in the scattering 
efficency. 
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I. INTRODUCTION 

In recent years, quantum dots (QDs) have emerged 
as a powerful tool to tailor the light-emission proper- 
ties of semiconductors i Applications range from quan- 
tum dot lasers and ultrafast amplifiers to cavity-quantum 
electrodynamics and nonclassical light generation.-^ As a 
new material system, group-Ill nitrides are of intense 
current interest due to their extended range of emis- 
sion frequencies from amber up to ultraviolet as well 
as their potential for high-power/high-temperature elec- 
tronic devices. '^•'^ While the nonradiative loss of carriers 
due to trapping at threading dislocations lowers the ef- 
ficiency of group-Ill nitride quantum- well light emitters, 
this effect is reduced by the three-dimensional confine- 
ment in QDs.^ Studies of the photoluminescence spectra?^ 
and dynamicsi*^ are used to demonstrate and analyze 
efficient recombination processes from the localized QD 
states. 

Extensive work has been done to study the electronic 
states in group-Ill nitridesi^ The valence-band structure 
has a strongly non-parabolic dispersion and a pronounced 
mass anisotropyjiS Nitride-based heterostructures with a 
wurtzite crystal structure are known to have strong built- 
in electrostatic fields due to spontaneous polarization 
and piezoelectric effects which have been analyzed in ab- 
initio electronic structure calculation o^^i^^ and in com- 
parison with photoluminescence experimentsii^ Tight- 
binding calculations of QD states have been used to study 
free-carrier optical transitionsi^ 

While the aforementioned theoretical investigations 
are devoted to the single-particle states and transitions, 
it is known from GaAs-based QD systems, that the emis- 
sion properties are strongly influenced by many-body ef- 
fects. Self-organized QD systems, grown in the Stranski- 
Krastanov mode, exhibit a single-particle energy spec- 
trum with discrete energies for localized states as well 
as a quasi-continuum of delocalized wetting layer (WL) 
states at higher energies. The carrier-carrier Coulomb 
interaction provides efficient scattering processes from 
the delocalized into the localized states (carrier capture) 



as well as fast transitions between localized states (car- 
rier relaxation), which can be assisted by carriers in 
bound (QD) or extended (WL) states. The dependence 
of the scattering efficiency on the excitation conditions 
has been calculated for GaAs-based QDs on various levels 
of refinement li^iiSiiLiSiiS The scattering rates are of cen- 
tral importance for the photoluminescence as well as for 
the laser efficiency and dynamics. The same interaction 
processes also lead to a renormalization of the electronic 
states (resulting in line shifts for the optical transitions) 
as well as to dephasing (line broadening) effects, which 
directly determine absorption and gain spectra:^*^ 

With the current attention on nitride-based QDs the 
question raises, to which extent previous results are mod- 
ified by the peculiarities of this material system. Specif- 
ically, we study how the strong electrostatic fields and 
the corresponding changes of the single-particle wave 
functions and energies influence the carrier-carrier scat- 
tering processes. For this purpose, we have to analyze 
the competing influence of the internal fields and of the 
many-body renormalizations. Since previous investiga- 
tions of carrier-carrier scattering in semiconductor QDs 
have been performed for free-carrier energies entering the 
scattering integrals, an independent second purpose of 
this paper is the inclusion of self-consistently renormal- 
ized energies. 

Renormalizations of the single-particle energies are due 
to direct electrostatic (Hartree) Coulomb interaction, ex- 
change interaction, and screening. While these effects 
generally contribute due to possible charging of the QDs, 
in the considered wurtzite structure they are addition- 
ally modified by the presence of the built-in fields. In 
this paper we show that the discussed changes of the 
single-particle energies have a much stronger impact on 
the carrier scattering processes than modifications of the 
single-particle wave functions. 

Many-body effects were also considered for the case of 
few excited carriers restricted to localized states of GaN- 
QD.^ In this regime, the emission properties refiect the 
multi-exciton states. In our paper, we are interested in 
kinetic processes for the opposite limit of high densities 
of carriers populating both QD and WL states (as typical 
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for QD lasers). 

For the evaluation of scattering processes, we use ki- 
netic equations for the carrier occupation probabilities 
which include direct and exchange Coulomb interac- 
tion under the influence of carrier screening as well as 
population effects (Pauli-blocking) of the involved elec- 
tronic states. Based on the single particle states and 
Coulomb interaction matrix elements, scattering pro- 
cesses are evaluated on the level of second-order Born 
appr oximat ion . ® 

The paper is organized as follows. In Sec. ^ we sum- 
marize the main ingredients of our theory which include 
the kinetic equations, the single-particle states and their 
renormalization as well as the interaction matrix ele- 
ments. Details regarding the WL states and screening 
effects are given in appendices. In Sec. lIIII the QD model 
is described and the numerical results for energy renor- 
malization and scattering times are presented. 



II. THEORY FOR CARRIER-CARRIER 
COULOMB SCATTERING 

To analyze the role of carrier-carrier scattering in 
nitride-based QD systems, we use a kinetic equation with 
Boltzmann scattering integrals^ The dynamics of the 
carrier population fu{t) in an arbitrary state v is de- 
termined by in-scattering processes, weighted with the 
nonoccupation of this state, and by out-scattering pro- 
cesses, weighted with the occupation according to 

= (1 - - /.^r • (1) 

The in-scattering rate is given by 
27r 

I'1,1'2,I'3 

X - fv2)fvj{eu - + £,.2 - e<^3)}> (2) 

where the first and second term in the square brack- 
ets correspond to direct and exchange Coulomb scatter- 
ing, respectively. The matrix elements of the screened 
Coulomb interaction, Wnv^v^v^-, are calculated in Sec- 
tion IIIBI The delta-function describes energy conser- 
vation in the Markov limit. Single-particle energies Ev of 
state V are renormalized by the Coulomb interaction as 
discussed in Sections III Dl and |ll El A similar expression 
for the out-scattering rate 5°"* is obtained by replacing 
/->!-/• 

Scattering processes described by Eqs. and (0) con- 
serve the total energy and separately the electron and 
hole numbers. As a result, the combined action of the 
discussed scattering processes will evolve the distribu- 
tion functions of electrons and holes towards Fermi-Dirac 
functions with common temperature where the QD and 
WL electrons (holes) will have the same chemical po- 
tential /ig (^/t). During such a time evolution towards 



quasi-equilibrium the relative importance of various scat- 
tering processes is expected to change via their depen- 
dence on the (non-equilibrium) carrier distribution func- 
tions for WL and QD states. A direct measure for the 
efficiency of various scattering processes can be given 
in the relaxation-time approximationiS where one cal- 
culates the characteristic time on which a small pertur- 
bation of the system from thermal equilibrium is oblit- 
erated. We use this method to investigate the influence 
of the built-in electrostatic field in nitride-based QDs on 
the carrier-density dependent scattering efficiency. 



A. Quantum-dot model system 

Recent progress in tight-binding and k • p models has 
been made in calculating QD electronic single-particle 
states including the confinement geometry, strain and 
built-in electrostatic field effects. The special features 
of the wurtzite QD structures have been addressed in 
Refs. 14,22,23'24'|, while zinc-blende QD structures have 
been studied in Refs. 2 3,24 25,26 . 27.28] . Our goals dif- 
fer from these investigations in the respect that, for given 
single-particle states and energies, we need to determine 
Coulomb interaction matrix elements in order to calcu- 
late many-body energy renormalizations and scattering 
processes. For this purpose we choose a simple represen- 
tation of wave functions for lens-shaped quantum dota^ 
which allows to separate the in-plane motion (with weak 
confinement for the states localized at the QD position 
and without confinement for the states delocalized over 
the WL plane) from the motion in growth direction with 
strong confinement. This leads to the ansatz 

$,(r) = ^f(^)e(z)u,(r) (3) 

where the WL extends in the plane described by g = 
{x,y). if and ^ are the envelope functions in this plane 
and in the perpendicular growth direction, respectively, 
and u are Bloch functions, v represents a set of quantum 
numbers with / for the in-plane component (including 
the spin), a for the z-direction, and b is the band index. 
In the following we consider an ensemble of randomly 
distributed identical QDs with non-overlapping localized 
states. The total number of QDs, N, leads in the large 
area limit to a constant QD density tt-qd = ViTaA^oo N/A. 

Regarding the dependence of the results on the choice 
of wave functions, we find that not so much the particu- 
lar form but rather the correct symmetry is of relevance. 
Even tough a more accurate treatment is expected to mix 
the in-plane and z-coordinates, the ansatz Q is preserv- 
ing the symmetry. 



B. Coulomb matrix elements 

The interaction matrix elements of the bare Coulomb 
potential v{r — r') = e^/(47reoe|r — r'|) with the back- 
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ground dielectric function e are given by 



K{r)KA^')v{r^r')^.Ar')^.dr). (4) 

This expression is further specified witli the help of 
Eq. (PI, the Fourier transform of the Coulomb poten- 
tial, and by introducing the in-plane Coulomb matrix 
elements with the two-dimensional momentum q, 



v: 



b,b' 



Limiting the calculations to the first bound state of the 
strong confinement problem (in z-direction), all the a 
indices above take only one value and will be dropped in 
what follows. Also, the band indices associated with vi 
and are the same as those of v and V2, respectively, 
so that only two of them have to be specified. Therefore 
Eq. igj reads 



V, 



b,b2 



1 



q 



(6) 



The obtained separation of integrals over in-plane and 
^-components greatly simplifies the computational ef- 
fort. For practical calculations we use solutions of a 
two-dimensional harmonic potential for the in-plane wave 
functions of the localized states and orthogonalized plane 
wave (OPW) solutions, discussed in Appendix 1X1 for the 
in-plane components of the delocalized states. Then the 
in-plane integrals in Eq. © can be determined to a large 
extent analytically for all possible combinations of QD 
and WL states. The calculation of the wave functions 
in growth direction, which is used to evaluate the in- 
plane Coulomb matrix elements, Eq. (jSj), is outlined in 
the next section. Screened Coulomb matrix elements arc 
obtained from Eqs. © and © according to the proce- 
dure described in detail in Ref. jl9| . 



C. Quantum-confined Stark effect 

To account for the built-in electrostatic fields of the 
wurtzite structure in the growth-direction, which gives 
rise to the quantum-confined Stark effect, we solve the 
one-dimensional Schrodinger equation^S*^ 



2mh dz' 



+ U\z) 



eiz) = E,iz)eiz), 



where the potential 



U''{z) = U',{z) + U'{z) + UU^), 



(7) 



(8) 



consists of the bare confinement potential in z-direction, 
Uq{z), as well as the intrinsic electrical and screening 
fields, Up{z) and Ug^^.{z), respectively. The latter is 
the electrostatic field due to the separation of electron 
and hole wave functions by the built-in field. Following 
Ref. [sJl we calculate the corresponding screening poten- 
tial from a solution of the Poisson equation for a set of 
uniformly charged sheets according to 



2eoe 



dz'[nz')\^-\eizT] k-^'i-(9) 



Equations 0-® have to be evaluated selfconsistently 
for a given total (QD plus WL) carrier density Ngys in 
the system. 



D. Hartree-Fock energy renormalization 

The Hartree-Fock (HF) contribution to the energy 
renormalization of an arbitrary state v (QD or WL) is 
given by 



= + A[^ 



HF 



(10) 



where is the free-carrier energy and the HF shift fol- 
lows from 



.HF 



(11) 



The first part corresponds to the Hartree (direct) term 
and the second part is the Fock (exchange) contribution. 
The equation is written in a general basis. For a system 
with local charge neutrality, like bulk semiconductors or 
quantum wells, the Hartree term vanishes while the Fock 
term leads to an energy reduction as used, e.g., in the 
semiconductor Bloch equationsi^ 

For the QD and OPW-WL states discussed in this 
paper, the absence of local charge neutrality leads to 
Hartree terms which are evaluated in Appendix ^ 
Regarding the quantum numbers, introduced in Sec- 
tioning we further specify the following notation: For 
the in-plane envelopes we use the two-dimensional mo- 
mentum k for the delocalized WL states and a = (m, R) 
for the localized QD states, where R is the QD position 
and the discrete quantum numbers for a particular QD 
are collected in m. The spin index is tacitly included in 
either m or k. 

Following Appendix^ we find that the Hartree shifts 
of the WL states vanish due to compensating contribu- 
tions from QD and WL carriers. 



AH - 
^&.k — 



0. 



(12) 



For a random distribution of QDs this is related to the 
spatial homogeneity restored on a global length scale and 
to the global charge neutrality of the system. For the 
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same reason, the Hartree shift of a locahzed QD state is 
only provided by states from the same QD while Hartree 
shifts due to carriers in other QD and WL states com- 
pensate each other, 

fc' ,m' 

From electrostatics one expects the same result, pro- 
vided that the WL is modeled as a constant area charge of 
opposite sign to the QD total charge. Then the constant 
part of the Fourier expansion of the Coulomb matrix el- 
ements in Eq. ©, i.e. the q = contribution, for the 
QDs balances the constant area charge from the WL. 

On the other hand, since probes the local charge 
density at the site of the QD due to the contributions 
of QD and WL carriers, one expects an influence of the 
WL on the QD Hartree energy shift. Locally on the QD 
length scale, the WL states are not homogeneous as a 
result of the QD presence. This causes a departure from 
the picture where the WL states contribute only in an 
averaged manner (via the q = term). Intuitively, one 
would expect that an increasing amount of carriers in the 
WL will start to screen the Coulomb interaction between 
the QD carriers. Following this picture we therefore re- 
place the bare Coulomb potential with the screened one 
in Eq. (jTS)), i.e., A^,„ /\f^^. In Appendices Eland O 
more support is given to this argumentation. 

In the exchange terms, the summation over the QD 
positions can be performed directly, since the associated 
QD phase factors disappear for the Coulomb matrix el- 
ements Vmm'mm' and V^kak as Seen from Eq. ljAl|l and 
Eq. (|A2|) . respectively. The resulting exchange energy 
shifts contain the QD and WL contributions, 

aF ^ _ yh,b fb 

'-^bjUi / J ^ mn\' mm' J rri' 

m' 

- E ^^k'^k'/k'. (14) 

k' 

^fc.k = ^ ^ ■ E ^km'km' /m' 

- E^ki^^kk'/^- (15) 

k' 

Since there is an area associated with the Coulomb ma- 
trix clement Vkm'km', the QD density ^qd — N /A enters 
in Eq. 

E. Screened exchange and Coulomb hole 

The HF Coulomb interaction provides only the first 
approximation for energy rcnormalizations; correlation 
contributions can lead to important corrections. A fre- 
quently used extension of the HF energy shifts for high- 
density plasma excitation is the screened-exchange and 
Coulomb-hole approximation.-^ On this level, the com- 
bined contributions of Coulomb exchange interaction and 



Coulomb correlations beyond HF to the energy rcnormal- 
izations are approximated with the screened exchange 
term (where the bare Coulomb potential in the Fock term 
is replaced by a screened one) plus an energy shift de- 
noted as Coulomb-hole contribution. The Coulomb hole 
self-energy reads^^ 

Yp'^{Ti,V2,hM) = ^<5(ri - r2)5{h - ^2) 

X [Ty(ri,r2,ti)-1/(ri-r2)], (16) 

with the statically screened Coulomb potential W . In a 
general eigenfunction basis, this leads to the Coulomb- 
hole energy shift 

A™(i) - ^E [^--'--'(^) - ^--'-A ■ (17) 

In Eqs. (|14() and H15() we therefore replace the bare 
Coulomb potential with the screened one and substitute 
the Fock energy shift with the screened exchange plus 
the Coulomb hole, A^^ + A^". In the hmit 

of low carrier densities, the screened Coulomb potential 
reduces to the bare one, the Coulomb hole vanishes, and 
we recover the HF result. 



III. RESULTS 

We start this Section by describing the model we use 
and the material parameters employed for calculating the 
one-particle states. These correspond to the unexcited 
system. The next step is to renormalize these states in 
a self-consistent way to include the influence of the car- 
rier population and the presence of the built-in field, as 
described in Section ^ These single-particle properties 
enter the scattering integrals, Eq. (0). 

For the discussion of the numerical results one should 
bear in mind the two main consequences of the built-in 
field. On the one hand, the self-consistent energies en- 
tering the Fermi functions and the energy-conservation 
are sensitive to the electron-hole separation induced by 
the built-in electrostatic field. On the other hand, the 
Coulomb matrix elements are changed due to modifica- 
tions in the wave-function shapes and overlapping. 



A. The model and its parameters 

In the following examples we consider an InCaN/GaN 
QD-WL system using typical InGaN parameters^ listed 
in Table ^ For the alloy, we have interpolated linearly 
the dielectric constant e, the isotropic electron mass rrie, 
the hole mass parameters Ai and the spin-orbit splitting 

Aso. 

A specific feature of the wurtzite structure nitrides 
is the strong mass anisotropy of the holes. The mass 
of the heavy hole (HH) in the z-direction, is given by 
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TABLE I: Material parameters used in the calculations. 



Parameter GaN InN Ino.2Gao.8N 



Eg[300 K] (eV) 


3.438" 


0.756" 


2.677 


AE, (eV) 






0.457 


AEh (eV) 






0.304 


e 


8.9' 


15.3'' 


10.2 


me (mo) 


0.2" 


0.07" 


0.174 


Ai 


-7.21" 


-8.21" 


-7.41 


A2 


-0.44" 


-0.68" 


-0.488 


A3 


6.68" 


7.57" 


6.858 


Ai 


-3.46" 


-5.23" 


-3.814 


As 


-3.40" 


-5.11" 


-3.742 


A,o (eV) 


0.017" 


0.005" 


0.0146 


Quantum well (nm) 






3.0 


Fqd-wl (MV/cm) 






1.5 


J'barricr (MV/cm) 






0.75 



"From Ref. M- 
''From Ref. IsJ . 



m* = Too/ 1^1 + ^3 1 with the free electron mass toq and 
the mass parameters Ai and ^3.^'^ For the in-plane mo- 
tion, a strong hybridization between the HH and the 
light-hole (LH) subband leads to nonparabolic bands. 
This is due to the fact that the small HH-LH splitting, 
induced by the spin-orbit interaction, is enhanced by nei- 
ther the strain nor the z-confinement. To include the hy- 
bridization effect we use a nonparabolic HH dispersion, 

= ~a{A2 + Ai)k' - ^Al+a^Alk\ (18) 

where a = /2mo and A2 = Aso/3. For the motion of 
the electrons, isotropic effective mass and parabolic dis- 
persion are considered. Conduction (AEe) and valence 
{AEh) band offsets for InGaN grown on GaN have been 
estimated by splitting the gap energy difference between 
bulk GaN and the alloy with a band offset ratio 60:40 for 
electrons and holes, according to Refs. 34-35]. 

Compared to usual zinc-blende structure, the wurtzite 
structure is characterized by large built-in electric fields. 
The strength of the fields is related to the sponta- 
neous polarization discontinuity at the heterojunction in- 
terfaces and the piezoelectric polarization. -'^^■-^2. Internal 
fields in InGaN/GaN heterojunctions of a few MV/cm 
have been reported,'^^'^^ typically in a sawtooth profile, 
where the field in the QD-WL system has a different 
magnitude and opposite direction to the field of the bar- 
rier (reflecting a set of capacitors with non-equal surface 
charges). For the field inside the QD-WL region and 
in the barrier, we consider i^QD-WL — 1-5 MV/cm and 
^barrier = 0.75 MV/cm, respectively. 

As the estimated effective hole massesiii are larger than 
the electron mass, we expect the QD's to confine more 
states for holes than for electrons. On the same basis. 



TABLE IL QD parameters used in the calculations. 



Parameter Electrons Holes 



Shells s, p s, p, d 

Level spacing (meV) 90.0 30.0 

£s (meV) -160.0 -80.0 

£p (meV) -70.0 -50.0 

£d (meV) -20.0 

QD density (cm~^) 10^° 



the level spacing which scales inversely with the mass, 
is larger for electrons than for holes. The in-plane QD 
confinement is modeled with a 2D parabolic potential, 
capable of binding two energy shells (s and p) for elec- 
trons and three (s,p,d) for holes. The degeneracies of 
these shells (apart from spin) are 1, 2 and 3 for s,p and d, 
respectively. For electrons we assume a level spacing of 
90 meV with the p-shell 70 meV below the WL contin- 
uum edge, while for holes we assume a level spacing of 30 
mcV with the d-shell 20 meV below the WL continuum 
edge. Thus the e., — fe.s QD transition is close to the range 
given m Ref. '8"37'|. The QD parameters are summarized 
in Tabic ^] Finally, we assume 2;-direction confinement 
wave functions which are band dependent but equal for 
QD and WL states. 

The harmonic oscillator (HO) inverse localization 
length (3 is deduced from the level spacing via Swho = 
h'^P'^/m*. For electrons the effective mass of Table is 
taken. For holes we use the mass resulting from Eq. H18|) 
in the small k limit, to^ = too/|^2 + Ai\. This is justi- 
fied by the typical QD diameters of lOOA - 200A , which 
correspond to the region around A: = where the first 
term of Eq. (|18|l is dominant. 



B. Schrodinger and Poisson equations 

The charge separation of carriers along the growth di- 
rection z under the influence of the built-in electrostatic 
fields, obtained from a self-consistent solution of Eqs. {T))- 
is shown in Fig.n(a) and (b). A clear overlap reduc- 
tion of the wave functions for electrons and holes is ob- 
served. This leads to a decreasing form factor (the double 
integral in Eq. (jSJ which modifies the 2D Fourier trans- 
form of the Coulomb interaction) for the electron-hole in- 
teraction while the stronger carrier localization increases 
the form factors for electron-electron and hole-hole in- 
teraction. Thus, the presence of the built-in field leads 
to an effective reduction of the electron-hole interaction, 
while the electron-electron and the hole-hole interaction 
is enhanced. The screening field turns out to be a small 
correction to the strong built-in field for the range of 
densities considered. 
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FIG. 1: Total confinement potential U{z) (solid line) and 
modulus amplitude of the z-component wave function 
(dashed line) along the growth direction z for electrons (a) 
and holes (b) in the presence of the built-in field at a total 
carrier density A'sys = 10^° cm~^. The dotted lines are the 
wave function amplitudes for zero field. 
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FIG. 2: Renormalized energies as a function of the total car- 
rier density in system A'sys for s-shell (solid lines), p-shell (dot- 
ted lines), d±-shell (dashed lines), and WL k = (dashed- 
dotted lines). Calculation with (a),(c) and without (b),(d) 
electrostatic field are shown for electrons (a),(b) and holes 
(c),(d). The temperature is 300 K. 



C. Renormalized energies 

For a total carrier density A^sys, the corresponding 
carrier distributions fi, in thermal equilibrium are used 
to determine the renormalized energies from the self- 
consistent solution of the equation Si, = + A^^-f A^^ + 
A^^, for the QD and WL states. For we use the en- 
ergies listed in Tabic HTl for the QD bound states and the 
dispersion law of Eq. (|18|) for the WL hole states. 

The density dependence of the renormalized electron 
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FIG. 3: Same as Fig. |5| for different contributions to the 
renormalized p-shell energy Sp which are the free energy Sp 
(dashed dotted line), screened Hartree shift A™ (dotted line), 
screened exchange shift Ap"'^ (solid line), and Coulomb hole 
shift A^" (dashed line). 



and hole energies is shown in Fig. 12 (a) and (c), respec- 
tively, in the presence of the built-in field. The corre- 
sponding results without the built-in field are given in 
Fig. 121 (b) and (d) for comparison. The renormalized d- 
shell is split into two degenerate d± states and a do state, 
with a separation of a few meV. 

Quantitatively, the QD levels experience a smaller en- 
ergy shift for the built-in field compared to the zero field 
case. The origin of this difference lies in the Hartree term, 
which reflects the electrostatic interaction of a given car- 
rier with all the others. The field-induced change of the 
^-confinement functions tends to separate the electrons 
from the holes and, as a consequence for both, the repul- 
sive part of the Hartree term is increased and the attrac- 
tive part is decreased. This effect is illustrated explicitly 
in Fig. O for the p-shell where the four different con- 
tributions to the renormalized energies are shown. For 
electrons the Hartree shift is repulsive both in the pres- 
ence and in the absence of the built-in field, but more 
so in the former case. For holes, the built-in field makes 
the Hartree term less attractive. In both cases the net 
result is a set of shallower bound states. For the sake of 
completeness we mention that the different sign of the 
Hartree field for electrons and holes comes from the dif- 
ference in the QD population of electrons and holes and 
from the band dependence of the Coulomb matrix ele- 
ments. 

The screened exchange and Coulomb hole terms are 
not significantly different for with and without built-in 
field. As the extended WL states are only renormalized 
by the screened exchange and Coulomb hole term, we 
find an overall negative energy shift, lowering the free 
spectrum by an almost k-independent shift (not shown) . 
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D. Capture and relaxation times 

To quantify the importance of the different scattering 
processes, we study their dependence on the total carrier 
density in thermal equilibrium. Using the relaxation time 
approximationiSi one can introduce a scattering time t^, 
for a given process according to 

= (19) 

which gives a characteristic time on which the system will 
return to its thermal equilibrium distribution if exposed 
to a small perturbation, i.e., = —{fv — Fi,)/tjj, with 
Fn being the thermal equilibrium distribution. 

The scattering times are changed by the built-in 
field through different competing mechanisms. On the 
one hand the matrix elements are modified (see Sec- 
tion lf f f I3|l , on the other hand the QD energies are pushed 
closer to the WL continuum fSection |fff Cp . 

For illustrative purposes, we first study the influence 
of the built-in field on the scattering times by using the 
free (unrenormalized) energies within the scattering inte- 
grals. This reveals the field effect solely on the Coulomb 
matrix elements via the wave-function changes. As an 
example we consider a capture process, where an elec- 
tron or hole from the WL is scattered into the QD while 
another WL carrier (electron or hole) is scattered to an 
energetically higher WL state, as well as the reverse pro- 
cess. (Both contribute to the scattering time according 
to the relaxation-time approximation.) Thus the outer 
index in Eq. (^) is a QD state while the three summa- 
tion indices in Eq. (|2Jl belong to the WL states for this 
example shown in Figure 0] 

First we discuss the density dependence. The capture 
time decreases with increasing carrier density, as more 
scattering partners become available. Furthermore, cap- 
ture times for electrons are slower than for holes, because 
the QD electron levels are placed energetically deeper 
below the WL continuum edge. For energy-conserving 
scattering processes, the excess energy of a WL electron 
which is captured to the QD must be transferred to an- 
other carrier from the WL in the present example. In 
this way a capture to an energetically deep lying QD 
state is associated with large momentum transfer for the 
WL carriers. As the matrix elements have a Gaussian de- 
pendence on the in-plane momentum, scattering to WL 
states with high momentum is suppressed. 

The changes in the capture times produced by the 
field in the case of unrenormalized energies are minimal. 
This proves that the competing trends described in Sec- 
tion IIII Bl are nearly compensating each other, with a 
slight dominance of the effect of electron-hole scattering 
reduction. 

For the calculation of the capture times in Fig. [3 based 
on WL assisted capture processes?^ the renormalized 
energies are included in the scattering integrals 5*"''°"*. 
Now the capture times become substantially shorter in 
the presence of the built-in field compared to the zero 
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FIG. 4: Capture times for the p-shell as a function of the total 
carrier density in the system A^'sys using the free (unrenormal- 
ized energies), with (solid lines) and without (dotted lines) 
built-in electrostatic field. 



Field No Field 




01 2345601 23456 
N3y3 [ 10"cm-2 ] N,y3 [ 10"cm-2 ] 

FIG. 5: Capture times as a function of the total carrier density 
in the system A'sys for s-shell (solid lines), p-shell (dotted 
lines), do-shell (dashed-dotted lines), and d±-shell (dashed 
lines) . 



field case. Thus, the energy separation of the QD lev- 
els from the WL edge plays now the dominant role in 
the scattering times. Even though for the built-in field, 
the capture processes are somewhat slowed down by the 
reduction of the electron-hole interaction, they are still 
faster compared to the zero field case, where due to the 
electron-hole interaction the QD levels are energetically 
deeper in the QD. 

Figure El shows the WL assisted QD relaxation times 
for processes where a QD electron (hole) scatters to a 
different QD electron (hole) state by means of a WL car- 
rier. Alternatively, a QD carrier performs a transition to 
the WL while another WL carrier scatters into a different 
QD stateii^ Thus the outer index in Eq. Q belongs to a 
QD state while two of the summation indices in Eq. 
correspond to WL states and one label to a QD state. 
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FIG. 6: Relaxation times as a function of the total carrier 
density in the system A^sys with labeling as in Fig. |^ 



traction in the Hartree terms, both effects working in the 
same direction of shallower confined levels. This in turn 
causes an enhancement of the scattering efficiency. 

As expected for the QD-WL system, the rates for car- 
rier capture and relaxation strongly depend on the den- 
sity of excited carriers in the localized and delocalized 
states. For intermediate densities, the scattering efh- 
ciency increases with carrier density. For large densities, 
Pauli-blocking and screening of the interaction matrix 
elements slow down a further increase of the scattering 
rates. For typical InGaN QD parameters, a QD density 
of 10^'^ cm~^ and a carrier density of lO^'^ cm~^ at room 
temperature, direct capture of electrons (holes) to ex- 
cited QD states results in scattering times on the order 
of 100 (10) ps. Relaxation times for scattering between 
the QD hole states are more than one order of magnitude 
shorter than capture, while at elevated densities the elec- 
tron p-shell relaxation is of the same order of magnitude 
as capture. 



Mixed QD relaxation processesjlSL where e.g. a QD elec- 
tron scatters out to the WL while another hole from the 
WL scatters down to the QD hole state play only a mi- 
nor role due to the charge separation of the electrons and 
holes caused by the built-in field. 

Generally, the relaxation times for holes are one or two 
orders of magnitude shorter than for electrons, since the 
QD energy level spacing is larger for the latter. With 
built-in field the relaxation times become shorter com- 
pared to the zero field case, where the QD energy level 
spacing is larger (with the exception of a slightly slower 
relaxation of the s-shell electrons). Relaxation times for 
scattering between QD states are in general more than 
one order of magnitude shorter than capture. For the p- 
shell electron relaxation, a saturation effect due to Pauli 
blocking is observed at higher densities, which leads to 
comparable capture and relaxation times. 

IV. CONCLUSION 

The presence of the built-in electrostatic field in 
wurtzite heterostructures causes a charge separation of 
electrons and holes along the growth direction, which in 
turn reduces the electron-hole interaction and increases 
the electron-electron and hole-hole interaction. 

Our results show that the usually discussed impor- 
tance of this effect on the interaction matrix elements 
only weakly influences the scattering rates. It turns out 
that the change of the self-consistently renormalized en- 
ergies due to charge separation leads to a much stronger 
modification of the scattering rates. Specifically, for the 
infiuence on the interaction matrix elements, the reduc- 
tion of electron-hole scattering is partly compensated by 
the increase of electron-electron and hole-hole scattering. 
In contrast, for the energy renormalization, the charge 
separation leads to increased repulsion and decreased at- 
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APPENDIX A: OPW STATES AND 
INTERACTION MATRIX ELEMENTS 

The basic idea of the following scheme is to construct 
an approximate single-particle basis for the combined 
QD-WL system which provides a feasible way to compute 
the interaction matrix elements, Eq. (0J. We start from 
localized QD states ipa , as introduced in Sections III Al 
and III Dl with a — (r7i,R), and WL states in the absence 
of QDs, which are assumed to have plane-wave envelope 
functions ^%i^Q) — in the WL plane with the 

two-dimensional carrier momentum k. Quantum num- 
bers for spin, band index, and confinement in z-direction 
are not explicitly written for notational simplicity. In 
the presence of the QDs the orthogonality condition of 
the basis is imposed by projecting the plane waves on 
the subspace orthogonal to the QD states, as outlined in 
Ref. . The WL states are therefore given by the OPW 
functions \ipw) = ^{I'pD-T^a \fa){'Pc\'pi))- Assuming 
QDs with nonoverfapping wave functions, the sum over 
a — (m, R) counts various QD states m at different QD 
positions R. For randomly distributed identical QDs, the 
normalization is given by — 1 — A^X^m l(¥''ral'/'k)P- 

This scheme allows to evaluate the in-plane integrals 
(i/|e'i-e|z/') = / d^g ipl{Q)e^'i-^ip^'iQ) which appear in 
the Coulomb matrix elements, Eq. ij^Jl, for various com- 
binations of QD and WL states. When i/ and i^' are two 
QD states, one obtains 

( a \e'^-^\ a') = {m \e'^-^\ m' ) e'^'^ 6r,r> , (Al) 
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with the QD positions R and R'. For combinations of 
QD and WL states, one finds 

( a |e"i-^| v'k' ) = ( m \e'^-^\ ) e*(k'+q)-R , (A2) 

and for two WL states 

( y^k W^-°\ ¥>k' ) -^.q+k' I?oPw(k,k',q), (A3) 

follows with 

1 



i:>OPw(k, k',q) 



TVkA^k' 



(A4) 



'■[1- nY.\{^'M? ~ nY,\{^M)? 



the law of large numbers, that the distribution of these 
quantities is sharply peaked around their configurational 
averaged value. Therefore we may replace the QD con- 
tribution in Eq. (|B1|I by 



A 



R,R' b' ,m' q 



X (m e 



|TO)(m'|e'i-e|m') 



X e 



-j:q-(R-R') 



(B3) 



Note that the resulting Hartree shift does not depend 
on the QD position any more. The summation over the 
random positions R, R' is evaluated as in the disordered 
system theory (see e.g. Reft2l): 



The orthogonality requirement of the wave func- 
tions is directly reflected in the interaction vertices 
{vi\e'^^'^'^\v2) = ^ui,i/2 foi' q = 0. Meaningful results 
for the interaction matrix elements can be expected only 
when the approximate model shows the same behavior. 
Since we start from orthogonal QD states at a given QD 
position and assume nonoverlapping wave functions for 
different QDs, Eq. I|A1() reduces to a Kronecker delta 
for q = 0. QD and OPW-WL states are orthogonal by 
construction, i.e., the requirement is also fulfilled for 
Eq. I)A2|I . As described in Ref. it is the assump- 
tion of randomly distributed QDs which, in the large- 
area limit restores - on average - translational invariance 
and provides mutually orthogonal OPW states such that 
the above requirement is also obeyed by Eq. IIA3I) . Note 
that i:)opw(k,k,0) = 1. 



^/(R)5(R')= E /(R)5(R')+E/W5W 

R^^^R' R 

= N^{f))-{{g))^-N{{fg)), (B4) 



R,R 



where {{Fjj — 1/A J cPRFCR) denotes the configuration 
average. The first term is the uncorrelated average of 
the two random variables, while the second takes into 
account that for the same point they are correlated. 
In our case /(R) = e-'i-^,g(R') ^ e^^'^', ((/)> = 
= (5q,o and ((/<?)) = 1, so that one may write: 



. H,QD 



E 



v: 



b,b' 



mm' m' m •> m 



(B5) 



APPENDIX B: HARTREE ENERGY 
RENORMALIZATION 

Starting from Eq. the Hartree energy shift of the 
QD states has contributions form the QD and from the 
WL carriers, 

'-^b,a — ^ aa'a'a J a' ^ ''ak'k'a Jk' 



fc',k' 



= A 



H,QD 



A 



H.WL 



(Bl) 



The QD contribution can be specified further using the 
notation introduced above and Eq. jnj, 

A^r = EE7E^^''(q)/- 



A 

R' h',m' q 

X (m|e-*i-s|TO)(m'|e'i-e|m') 

^ g-iq-(R-R')^ 



(B2) 



The result depends on the QD positions through the 
phase factor arising from the interaction vertices as given 
by Eqs. ljAip - (|A3p . In the large area limit one obtains, by 



The first term, arising from the Coulomb interaction 
between different QDs is proportional to the total QD 
charge density, while the second one describes the 
Hartree interaction inside a given dot. 

The WL contribution to the QD Hartree energy shift 
can be evaluated similarly. 



H,WL 



= EiE^'"'''(q)/- 

b',k' q 

X (TO|e-*i-^|m)e-*'i-^ 
X (5k',k'+qi^opw(k, k',q) 

= TE^''''(q = 0)/k'- 



(B6) 



b'.k' 



The result has the same structure as the first term of 
Eq. (|B5|I and adds the WL charge density to the QD 
contribution. By charge neutrality these terms cancel 
each other and, as expected, the Coulomb singularity at 
q = is removed. One is left with the second term of 
Eq. ||B5)| . which gives Eq. 

Following the same steps as above we find for the WL 
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Hartree energy shift 

A H _ yb,b' rb' , yb,b' rb' 



b'k' 



= V 



= 0, 



''^^'(q = 0){nQDE/™' + ^E/k'} 



h'm' 



(B7) 



which vanishes due to global charge neutrality. 





k' =k+q 



7!(R i;jR jTiR ;j!R 

(b) (c) 



APPENDIX C: WL SCREENING 
CONTRIBUTIONS TO THE QD HARTREE 
INTERACTION 

In our description of the Hartree terms in the previous 
Appendix ^ the summation over randomly distributed 
QDs restores, in the large area limit, the in-plane trans- 
lational invariance of the OPW-WL states. On this level, 
only the averaged QD properties enter - a picture which 
is consistent with the expectation that in a system with 
a macroscopic number of QDs, like a QD-laser, only the 
averaged properties of the QD ensemble should be impor- 
tant. On a local scale at a QD position, however, the WL 
states do not obey translational invariance since pertur- 
bations of the WL states due to the QD appear. In truly 
homogeneous systems one has a q = Coulomb singular- 
ity which is canceled out by the global charge neutrality 
and no other Hartree contribution is present. We have 
shown in Appendix ^ for the system of randomly dis- 
tributed QDs on the WL, that a similar singularity is 
produced by the configuration averaging and is removed 
by global neutrality arguments. Then only the localized 
QD carriers are subjected to Hartree fields induced by 
carriers in the same QD. 

In this appendix, we reexamine the result using the 
Green's function (GF) formalismi^ The more refined 
treatment shows that WL carriers can provide correc- 
tions to the QD Hartree shift which can be cast into the 
form of screening contributions. 

Some low order diagrams in the GF expansion, describ- 
ing terms of the Hartree contribution to the QD energies, 
are shown in Fig.[7| Since the QDs arc identical, the QD 
propagators are position independent, but the interac- 
tion vertices contain phase factors related to the position 
and to the adjoining momenta, as given by Eqs. (|A1|) - 
(jA3|) . The procedure described in the previous Appendix 
amounts to the averaging of the Hartree self-energies. 

To begin with, the self-energy in diagram (a) has the 
form 

R,R'6',m' q 

X (m|e-'i-e|TO)(m'|e*i-^|m') 

X e-^i (^-^') . (CI) 

In a self-consistent calculation, with the equal-time GF 
related to the population factors in the usual way, this 



FIG. 7: Hartree diagrams discussed in the text. 



leads to the QD Hartree term of Eq. ljB3p . In the large 
area limit, using Eq. (|B4p . one obtains the two terms of 
Eq. (jB5|) . As noted before, the first term corresponds to 
averaging the 'tadpole head' of diagram (a) as if it would 
be independent of the rest of the diagram, while the sec- 
ond term contains the contributions of the correlations. 
It is easy to see that diagram (b) of Fig. leads to the 
WL Hartree contribution spelled out in Eq. (|B6p . This 
term, together with the first term of diagram (a) contain 
the Coulomb q = singularity specific to homogeneous 
systems. The second term of diagram (a) is an example 
of a local field which is not averaged out, the source of 
this field being 'in phase' with the charges that probe it. 

A similar analysis is valid for diagram (c). Its phase 
factors are the same as in diagram (a) and again one 
has two contributions. The first one, coming from the 
uncorrelated averaging, leads to a self-energy containing 
only the q = contribution, 

X l/''-^'(q = 0)Gl,(i',t'), (C2) 

where summation and integration over the inner variables 
is assumed. Since this entails the restriction k = k', the 
diagram contributes to the renormalization of the WL 
propagator of index k in diagram (b). Therefore this 
is already included in a self-consistent calculation. More 
interesting is the second term, arising from the correlated 
averaging, 

^yt+im',ra:Ait'X)- (C3) 

In this case the summation over the momentum transfer 
q remains unrestricted. The structure is similar to the 
second term of diagram (a), i.e., it corresponds to the 
intra-QD Hartree field, but with the additional k and 
k' = k + q WL propagators forming a Lindhard loop. 
The loop describes the screening of the intra-QD Hartree 
field by the WL carriers. 

The usual procedure in the GF theory is to leave the 
Hartree interaction unscreened and to add the Lindhard 
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loop of diagram (c) to the 'tadpole head' GF of the di- 
agram (b). This avoids the double counting of such di- 
agrams. As a result, a nondiagonal (k, k') GF appears 
in the self-consistent Hartree loop of diagram (b). Al- 
ternatively, one can avoid double counting by keeping 
only momentum-diagonal WL propagators and leave the 
Lindhard loop for the screening of the Coulomb line. We 
have chosen this second approach, which also consider- 
ably simplifies the formalism. 

A fully systematic analysis of all the possible diagrams 



and the action of the configuration averaging over them 
is way beyond the scope of this paper. The approxi- 
mation proposed here includes the following physically 
important features. The random phases associated with 
the QD positions give rise to a q = singularity, which 
is canceled out by the global charge neutrality. On the 
other hand, the intra-QD fields are not influenced by the 
phase factors and therefore are not averaged out. The 
same is true for the local WL charges that respond to 
these fields and induce their screening. 
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